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Abstract 

We have computed the one-loop electroweak expression of diagonal and non diagonal stop- 
antistop and sbottom-antisbottom production from initial state gluons at LHC. We have inves- 
tigated the possibility that the one-loop effects exhibit a dependence on "extra" supersymmetric 
parameters different from the final squark masses. Our results, given for a choice of twelve SUSY 
benchmark points in the MSSM with mSUGRA symmetry breaking, show that in some cases a 
mild dependence might arise, at the percent relative level, of not simple experimental detection. 
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I. INTRODUCTION 



If Supersymmetry exists, and the super particles masses are not all unfairly large, LHC 
will be able to produce at least a fraction of these new creatures. In particular, the ex- 
pected values of the cross sections for squark-antisquark pairs production should allow their 
relatively quick identification. This might be particularly true for the case of final stops, 
that are supposed to be the lightest squarks in the available theoretical framework of the 
Minimal Supersymmetric Standard Model (MSSM). Not surprisingly, three calculations of 
stop-antistop production already exist, two at the electroweak Born level (including NLO 
QCD effects) for diagonal [1] and non diagonal [2] production and a very recent one [3] for 
diagonal production at electroweak NLO. For what concerns the dependence on the MSSM 
involved parameters, at the Born level for diagonal production this is limited to the masses 
of the two produced stops, conventionally defined as t\ (the lighter one) and t 2 . For non di- 
agonal production, the calculation of [2] is done for the electroweak s-channel Born diagram 
with Z exchange, whose value turns out to be (surprisingly) possibly larger than that of the 
(kinematically depressed) NLO QCD diagrams, and in principle potentially dependent on 
the stop mixing angle. Unfortunately, the predicted value of the cross section is in this case 
very small, and an experimental measurement does not seem to be easily performable, at 
least in a first LHC running period. In conclusion, the available calculations at the elec- 
troweak Born level for stop-antistop diagonal production only depend on the stop masses, 
while those for non diagonal production appear of non trivial experimental determination. 
Possible observable effects from supersymmetric parameters different from the stop masses 
might only arise at the next (NLO) electroweak one-loop order. A complete and exhaustive 
estimate of these NLO effects has been performed very recently [3], and several important 
features have been stressed. In particular, the calculation contains, beyond the one-loop 
corrections to the Born (LHC dominant) initial gluon-gluon state, the one-loop corrections 
to the Born initial quark- ant iquark state and also the contributions from the photon-induced 
gluon-photon fusion channel, with the inclusion of QED effects with soft and hard photon 
emission. Briefly, one discovers (a) that the one-loop corrections to the less Born relevant 
quark- ant iquark initial state can be, in some cases, competitive with those coming from the 
gluon-gluon state and (b) that the effect of the gluon-photon channel can be larger than 
those of NLO electroweak nature. These results have been derived for an illustrative set of 
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four benchmark points, labelled as SPSla, SPSla', SPS2 and SPS5, and are given for the 
integrated diagonal cross section. For what concerns the extra SUSY parameter dependence, 
the analysis of [3] has been performed for the SPSla' benchmark point, with the conclusion 
that it might only arise from the initial gluon-gluon state but, apart from singular (i.e. of 
threshold type) effects it would be numerically modest (at the percent level). 
The process of stop-antistop production is not the only third sfamily case considered in 
the literature. The production of sbottom-antisbottom, that might be an interesting source 
e.g. of very light stop subsequent decays, has also been studied in [2], together with the 
stop-antistop one, for both diagonal and non diagonal cases, also considering the mixed 
stop-sbottom production, at the Born level and for the two benchmark points SPSla and 
SPS5. In fact, in [2] one mSUGRA parameter (the scalar mass m or the fermion one m^j) 
is varied from its default benchmark value and the effect on the cross section is shown. In 
practice, varying these masses changes automatically the stop and sbottom masses, and the 
plotted changes of the various rates are a pure consequence of the latter variations. For what 
concerns the non diagonal cases, the considered production mechanism is Z/W exchange at 
Born level, and the relevant parameters are the stop and sbottom masses, with a possible 
extra dependence on the stop mixing angle that would deserve a deeper investigation. 
In conclusion, from the available literature one derives a picture of stop and sbottom pro- 
duction that can be summarized as follows: 

1. Diagonal stop-antistop production has been computed at complete NLO electroweak 
order in the MSSM for mSUGRA symmetry breaking [3]. A very mild dependence 
on SUSY parameters different from the stop masses has been found in the considered 
SPSla' benchmark point. 

2. Diagonal sbottom-antisbottom production and non diagonal stop and sbottom pro- 
duction have been computed at Born electroweak level [2]. In both cases the relevant 
parameters for the total rates are the final squark masses, possibly the stop mixing 
angle. 

Given these premises, this paper has two different purposes. The first one is that of search- 
ing possible cases where the negative conclusions of [3] might be evaded. With this aim, we 
have repeated the calculation of the NLO electroweak effects on the gluon initiated diagonal 
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stop-antistop production, extending the analysis to a larger set of benchmark points. Since 
the "extra" (i.e. different from the stop masses) NLO SUSY parameter dependence was 
excluded in [3] for the initial quark-antiquark state, we have not included it in our analysis 
in which we have limited our QED calculations to the derivation of the soft photon contri- 
bution. The second purpose is that of performing a NLO electroweak calculation, in search 
of extra parameter dependence, of some of the processes considered at Born level in [2], i.e. 
non diagonal stop-antistop and diagonal and non diagonal sbottom-antisbottom production 
(where the relevant masses are now the sbottom ones). Again, we have limited our analysis 
to the initial gluon-gluon state for the diagonal process and to purely soft photon QED 
effects. 

Technically speaking, the paper is organized as follows. In Sections 2 and 3 we give the 
necessary details of the NLO electroweak calculation, including the treatment of the soft 
photon contribution, trying to limit the presentation to the essential ingredients. Section 4 
contains the definition of the proposed observables and the various computed NLO effects 
for different choices of SUSY parameters. In Section 5, some tentative conclusions are finally 
presented. 

II. THE KINEMATICS OF THE PROCESSES gg -> t a t* b , b a b* b 

We discuss in details the case of stop pair production initiated by 2 gluons gg — > t a t* b . 
The kinematic of the process g g — > b a bl is completely analogous to the stop case, the only 
differences being the substitution of stop masses and mixing angles with sbottom ones. 

Physical (mixed) stops and antistops are denoted as t a , and t%, with a, b running over 1 
and 2. They are obtained from the chirality states U, i — 1, 2 standing for L,R as 

explicitly 

t\ = cos 6 t tL + sin 9 t tn t 2 = — sin d t t L + cos 6 t tR (2) 
The momenta, polarization vectors and helicities are defined by 

g(p g ,e(\ g )) +g (p> g ,e'(\' g )) - F a (p a ) +f*(p 6 ) (3) 

The gluon polarization vectors depend on the helicities as 
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(1) 



£(9 K 0; ^r-7f°) ^ = ( 0! ^-^°) (4) 

We use also the kinematical variables 

S = (Pg + p' g ) 2 = (Pa + Pb) 2 ,U= (Pg-Pb) 2 = (Pg-Paf , * = (Pg - Paf = (Pg ~ Pbf (5) 

with 

P ff = ^(i; 0,0,1) Pl = ^(i; 0,0,-1) (6) 

p a = (E a ;p sin 9,0, p cos 9) p b = (E b ;-p sin 9,0, -p cos 9) (7) 

_ s + ml-mj _ s + ml-mj _ r— — _ 2p 

E a - E b - ^= P-VE a -m a (8) 

The helicity amplitudes F\ gj \> g , computed from the Feynman diagrams listed in the next 
Section using the polarisation vectors of Eq. (4), will appear with various combinations of 
colours of the external particles. Firstly, one can write the colour structure in the form 



Fx g , K = {^ 9 ^fe(y)]+F A 2 9)A ,[^ + ^(|)] 

+F14^ + F t4^ + F5 t J ] U (9) 

where i,j running from 1 to 8 refer to the gluon colours and a, (3 running from 1 to 3 refer 
to stop and antistop colours. 



The polarized cross sections of the process g g — > t a tl (averaged over initial and summed 
over final colours) read 

da{\g,\' g ) _ p ^ 
dcos9 2048vrs^' X<3 ^ 1 ; 

col 

and the unpolarised cross section is 

d(T 1 da(\ g ,\' g ) 
dcos9 4 ^ dcos9 [ ' 
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The colour summation can be explicitly written as 



y \F Xq \' | 2 = 12|F A X y I 2 + -Wl y | 2 + -(|F A 3 y | 2 + |F A 4 y H 

/ J 1 A 9> A 9 I I A 9> A 9 ' Q 1 A 9> A 9 ' Q Vl A 9 > A 9 1 1 A 9> A 9 ' 7 

col(ija(3) 

-\-\l( y Fx gt \^Fx gt \^ "^Ag ,Ag -^Ag ,Ag ) g (-^Ag,Ag-^Ag,Ag ^AgjAg-^AgjAg) 

- -F A 3 A ,F 4 A ; +16F 2 A ,F A 5 A , + 8(F 3 A , F A 5 A , + F 4 A ,F A 5 A , ) 

Q A 9> A 9 A 9> A 9 A 9' A g A 9 > A 9 V A 9> A 9 A 9> A 9 A 9 > A 9 A 9> A 9 7 

+24|F A 5 9iA ,| 2 (12) 

The Born terms 

The Born terms exist only for "diagonal" stop-antistop pairs (a = b). They are given by 
4 diagrams shown in Fig (1): 

j^Born j^Born A _|_ j^Born A' _|_ j^Born B _|_ j^Born C (13) 

(A) s-channel gluon exchange: 

A Born a = [i/*^] (4^ (e.e^ (14) 



(A') 4-leg g l g j t a t* a diagram: 



A B J rn A ' = [\s i3 + d^}(4na s )(e.e')5 ab (15) 



(B) stop exchange in the t-channel: 



A B a r B = - ^ [^] {e. P ){e'. P )5 ab (16) 



(C) stop exchange in the u-channel: 



167TO% r V X i 



A B a r c = - [yy] (e-P)(e'.p)S ab (17) 

(we have used e'.p' = —e'.p and e.p' = —e.p). 

As one sees the Born terms only involve 2 invariant forms 

h = (e.p) (e'.p) I 2 = (e.e') (18) 
(and 4 colour components, C = 1,4), so that writing the invariant amplitude as 

A = N 1 (s,t,u)I 1 + N 2 (s,t,u)I 2 (19) 
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the helicity amplitudes are given by: 

F Xg ,y g = - ^XgX'gP 2 sin 2 6 N^u) + + A g A^) N 2 (s,t,u) (20) 
From Eqs. (10,12) and (20) one obtains the polarized Born cross sections 

da^jXii x ^ ™ 2 sP ( m i u 28 + 36 ^ 2 cos2 6 1 r2n 

dcos# 24s 1 s 2 jL (l-/3 2 cos 2 0) 2 J 1 ] 

da B °™(\ g ,-\ g ) _ na^ 28 + 36/? 2 cos 2 9 

dcos9 384s 1 (l-/3 2 cos 2 0) 2 J 1 1 

in agreement with the results of Ref.[2],[4]. Note that, at this Born level, a Born (++) = 

a Born ( ) and a Born (+-) = a Born (-+). 

It is useful, for later discussions of one loop effects, to emphasize the energy and angular 
dependences of the two types of polarized cross sections which are illustrated in Figs. (6,7). 

At low energy the dominant cross sections are the so-called, [5], Gauge Boson Helicity 

Violating (GBHV) ones er(++, ) of Eq. (21). This arises because the invariant form I 2 

has no threshold suppression factor, contrarily to Ii which vanishes like (3 2 near threshold. 
However at high energy the GBHV cross sections become mass suppressed like m~ Is 2 , as 
one can check from Eq. (21), in agreement with the general HC rule of Ref.[5]. Consequently, 
as one sees in Fig. (6), between threshold (2mj) and about 3mj , the stop pair is essentially 
produced through a(++, — ), whereas for higher energies (y/s > 3m^) it is dominated by 
the GBHC cross sections (<r(H — , — h) of Eq. (22). In Fig. (7) we have shown the corre- 
sponding angular distributions which appear to be also totally different in the two cases, 
larger for central angles in a(-\ — , — h), see Eq. (22), as opposed to forward and backward 

peaks in <r(++, ), Eq. (21). These various features will be essential for understanding 

the sensitivity to one loop effects in this process at LHC. 

As shown in Ref. [3], the stop pair can also be produced through the qq channel, Fig. (8), 
with a cross section 

-^- = ^sin 2 ^ (23) 
dcosO 18s v ; 

and through photon-induced mechanisms present only at NLO. Even if at LHC these pro- 
cesses are depressed as compared to the gluon-gluon one because of the smaller (quark) or 
non-existent (photon) PdF's, the authors of Ref. [3] have shown that the one-loop corrections 
can be numerically significant or even bigger than those of the gluon fusion initiated process, 
but essentially independent of extra (i.e. different from the stop mass) SUSY parameters, 
and for this reason we shall not include them in our analysis. 



III. ONE LOOP ELECTROWEAK CORRECTIONS TO gg ->■ t a t* b , b a b* b 
A. Stop pair production: g g — ► t a tl 

The one loop electroweak contributions come from counter terms (c.t.) and self-energy 
(s.e.) corrections to the Born terms, and from triangle and box diagrams. We use the 
on-shell scheme [6] writing first the c.t.+s.e. corrections as: 

For a = b: 

A Born + c.t. + s.e. A = j^orn A[1 + ^ ^ 
A Born + c.t. + s.e. A' = j^rn A' ^ + ^ (25) 



A Born + c.t.+s.e. B = A Born B^ + _ ^aa(t) j ^ 



A Born + c.t.+s.e. C = ^orn C ^ + _ (2?) 

and for a + b, using A^ rn A > A ' = A B hh ° rn A ' A ' 

A Born+c.t.+s.e. A = A Born AJ^ ^ 

A Born+c.t.+s.e. A' = ^Born A'J^ ^ 

A Born+c.t.+s.e. B ABorn Brr 7 ^ab(t) 1 xBorn Brry ^afc^Q l for\\ 

4* = A aa [5Z ba - 2{t _ ml) ] + A bb [5Z ba - 2{t _ ml) ] (30) 

tb t a 



ABorn+c.t.+s.e. C _ ABorn Cfrv ^a&W 1 ,j£orn Crrv ^afc( M ) 1 /qi \ 

A* " Au, L^ta - 2(M _ m 2)J + ^ L^ta " 2(M _ m 2)J ( 31 ) 

with the c.t. terms expressed in terms of stops self-energies 

5Z aa = -[ — ^-i] p2=m ^ (32) 
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and for a ^ b 



TL ha {ml ) 



5Z ba — — 2 %~ (33) 

m~ — ra- 
th t a 



SZ ba = -[5Z* ba + 5Z ab ] (34) 
the renormalized s.e. functions being given by 

± aa (p 2 ) = Z aa (p 2 ) - S aa (mf )-(p 2 - m\ )[%P] p2=m2 (35) 



and for a ^ b 



-.2 ™2 2 ™2 



p — m~ rr — m~ 
HUP 2 ) = ^UP 2 ) + — \^a{m\ ) + — \KM b ) (36) 

tb ta ta 

The needed S(p 2 ) functions are obtained from the various {qV), (qH), (qx) bubbles and 
from the gauge boson (V) and the 4- leg (SStt) tadpoles depicted in Fig. (2). 

Triangle and boxes corrections are shown in Figs. (3,4,5). They affect respectively 
each sector (A), (A'), (B) and (C) appearing in the Born case. In the s-channel one finds 
"left" and "right" triangles and in the t- and u- channels one has "up" and "down" ones. 
Contributions of sector (C) are obtained from those of sector (B) by symmetrization rules 
for the 2 gluons: interchange of momenta, polarization vectors and colours (p g ,e(X g ),i) 
and (p' g ,e'(X' g ),j). The 3 types of boxes can be identified through their (clockwise) internal 
contents (SSVS), (qqx<l) an d (SSHS) for sector (B), the above symmetrization rules 
giving the crossed sector (C); S refer to all possible scalar states. 

These electroweak corrections can also be classified into: 

- gauge terms due to internal exchanges of gauge bosons (V — 7, Z, W) and of charginos, 
neutralinos (through their gaugino components), 

Yukawa terms due to exchanges of Higgs bosons (if), and also charginos, neutralinos 
(now through their higgsino components). 

The contributions of these various diagrams to the helicity amplitudes are obtained after 
colour decomposition according to Eq. (9) and are expressed in terms of Passariono-Veltman 
(PV) functions. The numerical computation is then done with a dedicated C++ code ex- 
ploiting the LoopTools library [7] . 



A first check of the computation is obtained by observing the cancellation of the 
divergences appearing in counter terms, self-energies, triangles and boxes. For some 
parts these cancellations occur separately in each sector, but for other parts they involve 
contributions from several sectors as required by gauge invariance. 

Another type of check is provided by the high energy behaviour of the helicity amplitudes 
which has to satisfy a number of "asymptotic" rules. 

As already noticed in Sect. II, at high energy, neglecting masses the only surviving Born 
helicity amplitudes obtained from the addition of (A+A'+B+C) terms are the GBHC ones: 

f - a< = (4to<)( ^ )[ _^_ + _^_ ] (37) 



with 

% = l*» + d»\j) + j i3 = \f + d*'(y) - if*(j) (38) 

in agreement with the theorem given in [5], whereas the GBHV ones (with X g = X' ) are 
mass suppressed (vanish like m 2 /s). 

From the general logarithmic rules established in [8] , one expects the one loop virtual 
electroweak contributions to give, for final unmixed L, R states (before applying the mixing 
matrices R a i), the following corrections to the GBHC Born amplitudes: 

Fx g , -x B = F*°™[1 + (39) 



c «(l + 26c 2 w ) 8 _ 2 ^ _ a(m 2 t +ml) s 

~ 2 

^ = ^ 2ln W 2 ~ ln M^- ^Ml [ln M^ (41) 

m t ~ m b 

m t = -— m b = (42) 

sin p cos p 

in which one identifies the "gauge" and the "Yukawa" parts. M is a typical mass scale 
whose precise value does not matter at Log accuracy. 
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We have checked analytically (by taking the leading logarithmic expressions of the PV 
functions) that the various self-energy, triangle and box contributions reproduce the above 
expressions in both gauge and Yukawa sectors. 



We conclude this Section by briefly discussing the treatment of infrared singularities. As 
usual, QED radiation effects can be split into a soft part which is infrared (IR) singular and 
a hard part including the emission of photons with an energy which is not small compared 
to the process energy scale. In this paper, we have only included the soft part which is 
necessary in order to cancel any the IR singularities associated with the photonic virtual 
corrections. As we have outlined in the Introduction, since we are only searching for extra 
SUSY parameter dependence, we have not included the hard part of QED effects. 
We denote by A Born and A 1 loop any invariant helicity scattering amplitude evaluated at Born 
or one loop level. IR divergences are regulated by a small photon mass A. IR cancellation 
holds for every helicity channel separately and we checked it numerically by taking the A — > 
limit of our calculation. The real radiation factorizes on the Born amplitude leading to 



The universal correction factor 5 S takes into account the emission of soft real photons with 
energy from A up to A£'™ ax <C y/s [9]. In our analysis, we have fixed AE™ ax = 0.1 GeV. 

B. Sbottom pair production: gg ^b a b\ 

The treatment of the one loop corrections for the sbottom case is again analogous to that 
of the stop, but the particles involved in the loops are different, so the numerical results 
for the one loop contributions obtained in the stop production process cannot be trivially 
transposed to the sbottom case. In practice, all the expression given in the above section are 
to be "mirrored" substituting every top-tagged quantity with its bottom-tagged counterpart. 
For this reason we give only a brief overview of the main differences that arise between the 
two processes. 

Since the main parameters that controls the processes are the masses of the final state 
squarks, we start from some considerations about how such masses affect the observables 
we are going to analyze. As it is possible to see in Tab. (I), the masses of stop and sbottom 




<J a ) + 2^ Born A 



1 loop = j R finite _ 



(43) 
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squarks change within a wide range of values depending on the scenario considered, and 
the thresholds for the production of q a q% vary accordingly affecting the values of the cross 
section. Thus, since at tree-level the only difference is the sbottom masses instead of the 
stop masses in t — m 2 and u — m 2 , considering scenarios with not too different masses, the 
tre-level cross sections should be comparable. In any case, at high energy all the masses can 
be neglected, so the cross sections are identical to a great approximation. 

At one loop level, two type of differences appear: a) the different masses in the various 
propagators, b) the different couplings in gauge, SUSY gauge and Yukawa couplings. This 
second type can be very simply pointed out by comparing the Sudakov coefficients controlling 
the high energy behaviour: 

c c OU+26C&,) b _ 2 ^ _ a(m? + mg) s 

<V. - £rl*»W ~ ^ " £Sk [ln ^ ] (45) 

c «« = ^ {2ln w ~ ln2 ik ] ' (46) 

the only difference coming from the R part. 

Again this should give only a slight difference at high energy when mass effects are 
negligible. 



IV. RESULTS 



A. Stop pair production: gg —> t a tl 



Our starting observable for this process is the invariant mass distribution defined as 
da(pp^tJt + X) 



inv 



= J dxi dx 2 dcos9g(xi, //) g(x 2 , n) 

x ^"^f S(y^S - M inv ) , (47) 
a cost' 

where VS is the proton-proton cm. energy, M inv is the t a + tl invariant mass, 6 is the stop 
squark scattering angle in the partonic cm. frame, and g(xi, fj) are the distributions of the 
gluon inside the proton with a momentum fraction Xi at the scale \i. We have used the LO 
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PDF set CTEQ6L [10] with /x = wq a +mi b - As we already mentioned, we include soft QED 
real radiation in order to cancel IR singularities. For the 2-^2 + 7 (soft) process we can 
identify M inv with the partonic cm. energy y/s. The shift induced by hard QCD radiation 
has been previously estimated for tt production in [11] and found to be at the level of a 
few percents. Since our observables will be defined by integrating over a wide range of Mj nv 
values, such a shift will be irrelevant for our conclusions. 

For our purposes, we have considered the total rate a tot of the process defined by in- 
tegrating the distribution da/dM inv over the full range of invariant mass values, from the 
threshold +m,T, for the diagonal light squark production (tit*) and for the non-diagonal 
case (tit*, + t 2 t*). 

Our analysis has been performed for a choice of a large number of SUSY benchmark 
points. More specifically we have considered 12 mSUGRA inspired points: the eight SPS 
points (SPSla, SPSla', SPSla slope, SPS2-6) [12] which allow, as far as SPSla, SPSla', 
SPSla slope, SPS2 and SPS5 are concerned, a direct comparison with the results of [2], [3], 
the two SU1, SU6 ATLAS points [13] and two light SUSY scenarios LSI and LS2 discussed 
in [14]. In Tab. (I) we have listed the values of the chosen set of parameters: m , m 1( / 2 , A , 
tan (3 and sign \x. 

Our results are shown in the next Figures. We have tried to draw a limited number of 
curves, that contain all the information that seems more relevant to us. With this purpose, 
we have first shown in Figs. (9,10,11,12) the shape of the differential distribution da/dM- mv 
with the related relative effect for two representative points, chosen as LSI and SPS5, both 
for stop and sbottom production. It is possible to see that both for the stop and sbottom 
cases the relative effect is positive near the threshold, but drops to negative values in the 
high invariant mass region. The same feature persists in all the remaining considered points. 
This can be understood from the discussion of the various helicity amplitudes in Sec. (II). At 
large M inv , the helicity conserving amplitude dominates with its Sudakov negative correction, 
while at small M inv the helicity violating amplitude is the larger one and receives a positive 
correction in a narrow region near the production threshold. 

However, due to the different masses of stops and sbottoms and to the different particles 
involved in the loops, there are substantial differences between the two processes: in the 
stop case the positive relative effects in the very low mass region soon vanishes, approaching 
typically a -10% limit, while in the sbottom case it is possible to note that threshold effects 
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(the peaks and troughs in the low mass region) are more pronounced and produce a typically 
bigger positive contribution which drops slowly to different limits in the high mass region. 
As a consequence, in the stop case one may expect to find a rather small effect in the total 
rate due to the cancellation between the corrections in these two regimes; in the sbottom 
case, by contrast, it is not possible, a priori, to predict whether the total one-loop effect will 
be positive or negative and to what extent, therefore to analyze the corrections to the Born 
results the numerical evaluation is necessary. 

In Tab. (II) we show the numerical values of the total rates for the different benchmark 
points. To allow a comparison with other calculations, we also show the values of the lighter 
squark masses that are fixed by the SUSPECT [15] and FeynHiggs [16] codes that we used. 

Our search of extra SUSY parameter dependence has been performed in the following 
way. For each benchmark point, we have varied in turn one of the four conventional param- 
eters (tan/3, TO1/2, mo and A ) in a reasonable range, and computed the variable relative 
one loop effect and rate. For practical reasons we have only considered in the diagonal stop 
case the largely dominant t\t\ component, and have shown the value of the t\ mass which is 
generated by the variation of the chosen parameter. We anticipate, to shorten our presen- 
tation, that for diagonal stop production we shall only show in Figs. (13-18) the complete 
numerical results for those cases that seem to us reasonably meaningful, in particular that 
correspond to a total rate not below the 1 pb (extreme?) limit. This choice selects the 
set of LSI, LS2, SPSla, SPSla', SPSla slope and SPS5 benchmark points, but to perform 
a comparison with Ref.[3] we have also included the (perhaps academical) case of SPS2. 
For the remaining benchmark points, given the negligible value of their rates, we have only 
shown, for academic information, the dominant relative one loop effects in Fig. (19). 
Figs. (13, 14) show the results that we have obtained for the point SPS5, which is perhaps the 
most relevant one. As a general feature, common to all the considered cases, one sees that 
the SUSY one loop effects are almost systematically negative and small, of the few percent 
size. For what concerns the dependence on the chosen parameter, one sees for SPS5 that 
the variation of m\ii can produce a maximal variation of the relative effect of approximately 
three percent. A smaller variation, of approximately 1.5 percent, is generated in the con- 
sidered range of tan (5. Varying m and A has essentially no practical effect (~ below one 
percent) on the one loop contribution. The latter remains, in all cases, of the few percent 
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at most. 

In Figs. (15-18), to save space, we have only shown the maximal relative variation and the 
corresponding parameter. This choice selects tan/3 for SPSla, TOi/ 2 for SPSla', tan/3 for 
SPS2, m 1 /2 for LSI and LS2. In the SPSla' case we have also plotted the tan/3 dependence 
to perform a comparison with Ref. [3]. In all cases, the relative one loop effect is negative 
and small, typically of the one-two percent size. 

A special case is that of the benchmark point SPSla slope, where the parameters m^, mo 
and A are related. In this case, we have plotted in Fig. (18) the variations with m 1 / 2 and 
tan (3. One sees that in the first case the relative negative effect can vary between one and 
four percent, remaining often in the three-four percent range. This represents the most rele- 
vant extra parameter dependence of our stop analysis. Varying tan /3 can produce a smaller 
(~ 1.5) effect, with an overall negative relative contribution in the five percent region which 
a priori might be visible with a dedicated experimental search. 

An important step is now the comparison with previous results. Concerning the total rates, 
one can see from our curves that the values obtained for the points SPSla, SPSla', SPS2 
and SPS5 essentially reproduce, taking the corresponding stop mass values, the gluon-gluon 
component of Ref. [3] Table 1. For the parameter dependence, Ref. [3] shows the SPSla' 
case but uses, apart from tan /3, a different set of parameters. A comparison of the tan (3 
dependences for this point shows a qualitative agreement, i.e. a small and negative effect 
that increases with tan /3, although our values are slightly larger, in the three percent range. 
We conclude in this case, in full agreement with Ref. [3], that the dependence on the extra 
SUSY parameters is for SPSla' extremely small. 

Another comparison can be performed for the SPSla slope and SPS5 cases with the plots of 
Ref. [2]. Again one can see an essential agreement between our one- lop results and the Born 
results of Ref. [2], as one would expect given the smallness of our one loop effects. 
The conclusion from our analysis of diagonal stop antistop production is that supersym- 
metric contributions due to extra SUSY parameters exist, but are generally apparently too 
small, at the few percent level, to produce an appreciable effect under realistic LHC ex- 
perimental conditions, at least in a first luminosity phase. Our next step has been that of 
repeating our analysis for the diagonal sbottom-antisbottom production. Here we have only 
considered the LSI and LS2 points, that would have a rate of the pb size. The results of 
our calculation are shown in the next Figures, that we now briefly comment. 
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As one sees from Figs. (20-23), the dependence of the effects on m and A is essentially 
negligible. For m 1 / 2 there is also no dependence on LSI, and a larger but irregular depen- 
dence (same values for different mi/ 2 ) on LS2. The dependence on tan f3 exhibits a different, 
and possibly appreciable, feature. One sees that the negative effect regularly increases with 
tan /3, like in the stop cases, but changing more, i.e. from ~2 percent to ~6 percent in the 
explored range. In particular, we believe that a relative effect of approximately six percent, 
in correspondence to a rate of the 5 pb size, might be, in principle, proposed for a highly 
dedicated experimental search. 

In conclusion, a reasonable picture that seems to emerge from our combined analysis 
of the stop-antistop and sbottom-antisbottom diagonal production processes is that of a 
possible, although mild, dependence of the one-loop electroweak effect in the mSUGRA 
scenario on extra parameters. Keeping this result in mind, we have also computed, for 
all benchmark points, the non diagonal total rate derived from one-loop gg electroweak 
diagrams. We remind the reader that a calculation of the non diagonal rate, derived from 
qq annihilation at Born level via Z exchange, already exists [2] for SPS5. In the stop case, 
the value that is obtained is larger than that coming from the kinematically depressed NLO 
QCD diagrams, and is equal to ~ 6 • 1CT 4 pb. In the sbottom case, the value that is obtained 
is equal to ~ 1.5 • 1CT 5 pb. In Tabs. (111,1V), we show the values that we have derived for 
the different benchmark points both for the one-loop gg diagrams and from the Z exchange 
calculations. One sees that the one-loop electroweak values are of the same size as those due 
to Z exchange and in some cases larger. This could have some relevance for the meaningful 
benchmark points. For example, in the LS2 stop production case, summing the one-loop 
with the Z-exchange contribution, one would get a total rate of approximately 10~ 2 pb. 
This is a factor 15 larger than the SPS5 point of [2], but realistically hard for experimental 
detection. A similar conclusion might be drawn for the rates of the remaining meaningful 
points if one sums the one-loop with the Z-exchange contributions. The results we have 
obtained for the sbottoms are similar, but because of the tiny cross sections involved, the 
experimental confirmation of our predictions will be again problematic. 
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V. CONCLUSIONS 



We have devoted our analysis to the search for extra (i.e. different from the final squark 
masses) parameters dependence in the processes of diagonal and non-diagonal stop-antistop 
and sbottom-antisbottom production from the gg initiated channel at EW NLO at LHC. 
With this aim, we have chosen twelve representative mSUGRA benchmark points and per- 
formed a variation of the mSUGRA parameters. We have verified in all cases the presence 
of a small (at the few percent level) relative difference of the effects with a more important 
role apparently played by different parameters for different points, in particular by tan j3 in 
a case of sbottom production. 

Certainly, the possibility of experimental verification of our conclusions would require 
very high luminosity scenarios and accuracies, representing a real challenge for the LHC 
experimental groups. This might, though, become interesting in case of a previous LHC 
super symmetric production, which might justify the idea of the dedicated experimental 
effort that we have mentioned. In this respect, we should mention that the possibility of a 
determination of SUSY parameters dependence from the process of stop-chargino production 
has been already considered by us at Born level in a previous paper [17]. In view of the 
results obtained in this present search, we are now considering the derivation of the EW 
one-loop effects on stop-chargino production. The calculation is already in progress. 
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TABLE I: mSUGRA benchmark points and masses of stops and sbottoms (all the values are in 
GeV) 
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TABLE II: Total cross-sections (in pb) for diagonal stop and sbottom production. The point SPSla 
slope has not been included since it coincides with SPSla at = 250 GeV. 
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TABLE III: Total cross-sections (in pb) for non-diagonal stop production starting from qq and gg. 
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TABLE IV: Total cross-sections (in pb) for non-diagonal sbottom production starting from qq and 
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Tree diagrams 




FIG. 1: Tree level diagrams for diagonal production g g — > t\ t\. 
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s.e. diagrams 




FIG. 2: Self-energy (generic) diagrams for diagonal production g g — > t\t\. They are composed of: 
scalar and vector tadpoles where the particles can be higgs bosons, sleptons and squarks (1) and 
SU(2)xU(l) gauge bosons (2); scalar, fermion and scalar-vector bubbles where the particles can 
be quark- x (3), squark-higgs (4) and squark-e.w. gauge boson (5). 
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up /down/left /right triangles 




FIG. 3: Up, down, left and right triangle (generic) diagrams for diagonal production gg — ► tit*. 
In the s-channel diagrams (1), (2), (3) we have labelled the internal gluon as a generic vector, while 
all the other vector particles are intended to be SU(2)xU(l) gauge bosons. Fermion loops and 
scalar loops involve quarks-x and squarks-higgs respectively, with the exception of diagrams (16), 
(17), (18) and (19) where the loops involve only quarks and squarks. 
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4 legs triangles 




FIG. 4: Four legs triangle (generic) diagrams for diagonal production g g — ► t\t\. As in the 
previous figure we label s-channel internal gluons in diagrams (1) and (2) as vectors, while all the 
other vectors are e.w. gauge bosons. The s-channel scalar in diagram (11) can be any neutral 
higgs. 
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Boxes 




FIG. 5: Box (generic) diagrams for diagonal production gg — > t\t\. Every vector is an e.w. gauge 
boson; 4-fermions boxes are made of 3 quarks and a %; 4-scalars boxes are made of 3 squarks and 
a higgs boson. 
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FIG. 8: Born diagram for non-diagonal squark production qq — > q a ql, (a / 6), via Z boson 
exchange. 
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FIG. 10: LSI, Born and one-loop distribution da/dM- mv for sbottom production. The right panel 
shows the percentual relative effect. 
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FIG. 12: SPS5, Born and one-loop distribution da/dM- mv for sbottom production. The right panel 
shows the percentual relative effect. 
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SPS5: scan over tan (3 
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FIG. 13: SPS5: scans over the mSUGRA parameters tan/3 and for diagonal stop production. 

The top panels show the percentual effect on the integrated cross section, the bottom panels show 

the variation in the value of the total cross section; the numbers above the curves in the bottom 

panels represent the value of the stop mass mj (in GeV) . 

34 




100 110 120 130 140 150 160 170 180 190 200 

m [GeVj 



SPS5: scan over A Q 

stop production 

-i | i | i | i | i | i | i r 




3 l l i l i l i i i l i l i l i l i l i l i l 

-1000 -990 -980 -970 -960 -950 -940 -930 -920 -910 -900 

A„[GeV] 



FIG. 14: SPS5: scans over the mSUGRA parameters tuq and Aq for diagonal stop production. 

The top panels show the percentual effect on the integrated cross section, the bottom panels show 

the variation in the value of the total cross section; the numbers above the curves in the bottom 

panels represent the value of the stop mass (in GeV) . 
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SPSla: scan over tan (3 
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FIG. 15: SPSla and SPS2: scan over the mSUGRA parameter tan/3 for diagonal stop production. 

The top panels show the percentual effect on the integrated cross section, the bottom panels show 

the variation in the value of the total cross section; the numbers above the curves in the bottom 

panels represent the value of the stop mass (in GeV) . 
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LS 1 : scan over m 
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FIG. 16: LSI and LS2: scan over the mSUGRA parameter for diagonal stop production. The 

top panels show the percentual effect on the integrated cross section, the bottom panels show the 

variation in the value of the total cross section; the numbers above the curves in the bottom panels 

represent the value of the stop mass (in GeV) . 
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FIG. 17: SPSla': scan over the mSUGRA parameters tan/3 and for diagonal stop production. 

The top panels show the percentual effect on the integrated cross section, the bottom panels show 

the variation in the value of the total cross section; the numbers above the curves in the bottom 

panels represent the value of the stop mass m^in GeV). 
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SPSla slope: scan over tan (3 
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FIG. 18: SPSla slope: scan over the mSUGRA parameter tan/3 and m^ 2 for diagonal stop 

production. The top panels show the percentual effect on the integrated cross section, the bottom 

panels show the variation in the value of the total cross section; the numbers above the curves in 

the bottom panels represent the value of the stop mass (in GeV). 
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FIG. 19: Dominant parameter dependence on one loop effects for the benchmark points with small 



cross section for diagonal stop production. 
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FIG. 20: LSI and LS2: scan over the mSUGRA parameter mo for diagonal sbottom production. 

The top panels show the percentual effect on the integrated cross section, the bottom panels show 

the variation in the value of the total cross section; the numbers above the curves in the bottom 

panels represent the value of the sbottom mass (in GeV). 
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FIG. 21: LSI and LS2: scan over the mSUGRA parameter Aq for diagonal sbottom production. 

The top panels show the percentual effect on the integrated cross section, the bottom panels show 

the variation in the value of the total cross section; the numbers above the curves in the bottom 

panels represent the value of the sbottom mass m^(in GeV). 
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FIG. 22: LSI and LS2: scan over the mSUGRA parameter mw2 for diagonal sbottom production. 

The top panels show the percentual effect on the integrated cross section, the bottom panels show 

the variation in the value of the total cross section; the numbers above the curves in the bottom 

panels represent the value of the sbottom mass m^(in GeV). 
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LSI and LS2: scan over tan (3 
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FIG. 23: LSI and LS2: scan over the mSUGRA parameter tan (3 for diagonal sbottom production. 
The top panels show the percentual effect on the integrated cross section, the bottom panels show 
the variation in the value of the total cross section; the numbers above the curves in the bottom 
panels represent the value of the sbottom mass m~ h (in GeV). 
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